Dynamic approach to fusion of massive nuclei. 
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Abstract 



The role of the entrance channel in the fusion-fission reactions leading to 
nearly the same superheavy compound nucleus is studied in the framework of 
dynamic model. The calculations are done for ^'^Ca + ■^'^^Pu and ''^'''^Ge + 
208p^ reactions which could lead to formation of superheavy element Z = 114. 
It is shown that for these reactions there is an energy window for the values 
of the bombarding energy at which a capture probability is sufficiently large. 
Together with the restriction coming from the intrinsic barrier for fusion, it 
helps to find an optimal value of the bombarding energy for a given projectile- 
target combination. 
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I. INTRODUCTION 



The cross section for the production of the superheavy elements depends on the choice 
of the projectile-target combination and the bombarding energy -E'c.m.- The optimal choice 
is determined by the requirements to have a larger fusion cross section and larger survival 
probability of a compound nucleus relative to fission. For a given projectile-target combi- 
nation, a larger value of the bombarding energy is needed to overcome the reaction barrier 
which is determined by the nucleus-nucleus potential and the dynamic barriers if they exist. 
However, the excitation energy of the compound nucleus increases with the bombarding 
energy. It decreases the survival probability relative to fission of a nucleus produced in a 
reaction and therefore puts a restriction on the upper value of a bombarding energy. To 
determine the optimal value of Ec.m. it is necessary to analyse a dependence of a partial 
fusion cross section which is proportional to a capture probability, on a bombarding energy. 
To do it, we require in dynamic model to describe the initial stage of a heavy ion collision. 
Such a model has been developed in our earlier papers [1,2] and it is the aim of the present 
paper to apply this model to calculate the capture probability. The latter quantity is deter- 
mined by the dynamic aspects of the reaction mechanism and by the depth of the pocket in 
the nucleus-nucleus interaction potential. We have calculated a nucleus-nucleus interaction 
potential using a double-folding procedure with the Migdal's effective forces [3]. As the 
examples we consider below the following reactions: ^^Ca -|- ^"^^Pu and ^"^'^^Ge -|- ^°®Pb. 



II. BASIC FORMALISM 

The cross section of production of the evaporation residues ((Ter-) 

oo 

aer{E) = Y.{21 + l)a{''%E, l)Wsur{E, I) (1) 
z=o 

is determined by the partial fusion cross section {ai^^{E)), 

a(^%E) = ar'''''%E)PcNiE, I), (2) 

capture / rri\ ^capture / rri\ /o\ 



where A is a wavelength, Pcn{E,1) is a factor taking into account a decrease of the fusion 
probabihty due to dinuclear system break up before fusion, -ppp^^'^^j^^^ jg capture prob- 
abihty which depends on the coUision dynamics and determines the amount of partial waves 
leading to capture. The cross section of production of the evaporation residues aer depends 
as well as on the probability {Wsur{.E,l)) that the compound nucleus survives during the 
deexcitation cascade at the bombarding energy E. 

To calculate capture probability -p^P*"^*^)^^) shall use a dynamic approach developed 
in [1,2]. In this approach, the system of equations is derived to describe the radial motion 
of colliding nuclei and an evolution of their intrinsic states during the heavy ion collision. 
The relative motion coordinate R(t) and the velocity R (t) are determined by solving the 
equations of motion 

^^{n{t)) R, +Y.ikAm] it) = -^^^^^ (4) 

where 

/i(R) = mATAp/{AT + Ap) + 5/i(R), 

5/i(R) is the dynamic contribution to the reduced mass /i, 7fcj[R(t)] is the friction tensor, 
iy(R) = V^(R) + 6V(Il) is the nucleus-nucleus interaction potential and 6V(Ii) is the dy- 
namic contribution to a nucleus-nucleus potential which is due to the rearrangement of the 
densities of the interacting nuclei during reaction. To calculate 5/i(R), 7fcj[R(t)] and 6V(Ti), 
it is necessary to find the occupation numbers of the single particle states. Since the excita- 
tion energy of the interacting nuclei changes significantly during the course of the collision, 
it is necessary to take into account the time dependence of the occupation numbers. An 
evolution of the occupation numbers has been defined by a numerical solution of the von 
Neumann equation for the single particle density matrix h with the Hamiltonian H which 
takes the following form in the second quantized representation 

i^(R(t),0 = Y^^p^p^p + IZ^Ta^a^ + Y,Vii,{'R{t))atai, + Ke.- (5) 

P T 

Here ^ is the short notation for the relevant intrinsic variables, the third term on the r.h.s. 
of the Eq. (^) can be written as 
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Viii'R{t))atai, = AS'(RW)apap, + Y ATT'(R(^))aJa2., + (6) 

P,P' T,T' 

EgprW^))(4aT + h.c.), 

where P = {np, jp,lp,mp) and T = {ut, jr^Ty^r) are the sets of quantum numbers 
characterizing the single particle states in the noninteracting projectile and the target nuclei, 
respectively. The last term on the r.h.s. of Eq. (|^) represents the residual interaction. Since 
an explicit allowance for the residual interaction is very complicated it is customary to 
take into account a two-particle collision integral in the linearized form (r-approximation) 
[1,2,4-6] 

= [n{K{t))Mt)] - '-^im - r^(R(t))] , (7) 

where h'^'^{'R{t)) is the local quasi-equilibrium distribution, i.e. a Fermi distribution with 
the temperature T{t) corresponding to the excitation energy at the internuclear distance 
R(t). All formulae needed to calculate 7fcj[R(t)] and 6V(li)) are given in [2,4,6]. 

The nuclear part of a nucleus- nucleus potential V"(R(t)) is calculated using the double- 
folding procedure between the effective nucleon-nucleon forces feff[p{x)] suggested by Migdal 
[3] and the densities of the interacting nuclei taken in the Woods-Saxon form 

|r - R^-(t)| - i?oi^(l + /3f Vaol^A')) 



pS^)(r,R;,(t),0;,,#)) 



1 + exp(^- 



(8) 



where Hk are the center of mass coordinates and Rqk are the half density radii of interacting 
nuclei K = 1,2; I32{K) are the quadrupole deformation parameters determined by the B{E2) 
to the first-excited 2+ state (its value is taken from [7]) and 6k are the axial symmetry axes 
orientations relative to R(t). Thus, we have a possibility to consider fusion at different 
mutual orientations of the interacting nuclei. 

The competition between complete fusion and quasifission of a dinuclear system formed 
after capture and its further evolution are described using the method developed in [8]. 
This method is based on the assumption that dinuclear system formed in the collision of 
two nuclei evolve to fusion by increasing its mass asymmetry. It means that the mass 



asymmetry degree of freedom rj = {At — Ap)/ {At + Ap) is the main dynamic variable. The 
internuclear distance R{t) takes the value corresponding to the location of the minimum of 
the nucleus-nucleus interaction potential for every given value of rj. The evolution of the 
system along mass asymmetry degree of freedom is described by a driving potential U {Z, I) 
which is calculated as 

U{Z, l) = B, + B2 + Uu{Rm) - Bo. (9) 

Here, Bi and B2 are the binding energies of the nuclei in a dinuclear system, Ui2{Rm) is 
the value of the nucleus-nucleus interaction potential at the minimum, Bq is the binding 
energy of the compound nucleus (the binding energies B^ are obtained from [9] and from 
[10] particularly for the superheavy elements). Therefore, a dinuclear system to be fused 
should overcome the intrinsic barrier {Bj^^) which is determined by the difference between 
the values of a driving potential located at the Businaro-Gallone point {rj = tibg) and the 
initial point corresponded to reaction under consideration. For the reactions considered 
below, the initial value of rj is smaller than tibg- The quasifission, which is in competition 
with the fusion is considered as a motion in the nucleus-nucleus interaction potential W{R). 
Thus, for quasifission, it is necessary to overcome a barrier of W{R). The competition of 
fusion and quasifission is taken into account by the factor Pcn{E,1), which is calculated 
using the following relation derived from the statistical model arguments 

p p(^* - 

'''' PiE* - B}^,) + p{E* - B^fY ^ ^ 

Here p{E* — B^) is the level density 



_ s.,) = ^(^-^X'" e.p^2.V,(e.)(^--B-)/6j ^^^^ 



2\M(eF)^2(eF) 



lg{eF){E* - B*)\' {E* - B*j,)V^ 

In Eq. (|10|), Bj^g is the barrier of the driving potential U{Z,l), which should be overcome 
on the way from the initial value of 1] to rj = 1. The B*j is the barrier of the nucleus- nucleus 
interaction potential which should be overcome if dinuclear system decays in two fragments, 
E* is an excitation energy of the compound nucleus which is equal to difference between 



Ec.m. and the mmimum of nucleus-nucleus potential {E* — Ec.m. — U{Rm)), gi,2{^F) are 
the single particle level densities of the fragments of the dinuclear system and g{e) is their 
sum, Krot is a factor taking into account rotation of a dinuclear system 

K...^^^^El^j,, (12, 

TT 

where J± is the rigid body moment of inertia for rotation around the axis perpendicular to 
the line connecting the centers of fragments. 



III. RESULTS AND DISCUSSION 

We consider below the following reactions which are discussed now as possible ways to 
search for superheavy element with Z = 114. They are ^^Ca + ^^^Pu (I) and ^^'^^Ge + ^^^Pb 
(II,III). 

Basing on the dynamic model developed in [1] (which is described concisely in the preced- 
ing section) we have calculated the capture cross section a^^*^^^(E) for the reactions under 
consideration. The results are shown in Fig. 1. It is seen that, for these reactions there is an 
energy window for the values of the bombarding energy at which a capture cross section is 
large enough to have a physical interest. The lower limit for the bombarding energy {Emin) 
is defined by a total nucleus-nucleus interaction potential W{R) = V{R) + 5V{R). Note 
that E^in is somewhat larger than the value of the entrance Coulomb barrier, because of 
the kinetic energy loss due to friction. So, £'mm is determined by a dynamic calculation. 
The upper limit (Emax) comes from an incomplete dissipation of the relative kinetic energy. 
Thus, the values of E^^n and E^ax are determined by the depth of the pocket in the poten- 
tial W{R) (Fig. 2) and by dissipative forces. If a bombarding energy is larger than E^nax 
the dissipative forces could not provide a complete dissipation of the relative kinetic energy 
and dinuclear system decays into two fragments instead of being fused. As it is seen from 
Fig. 1, reaction with the hghter projectile (I) has a larger value of the capture cross section 
than other two reactions (II) and (III). The reason is that for ^^Ca -|- ^^^Pu reaction the 
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pocket of the nucleus-nucleus interaction potential is deeper and wider than for ' "^Ge + 
208p^ (see Fig. 2). The potentials presented in Fig. 2 are calculated taking into account a 
deformation of the interacting nuclei assuming the tip-tip orientation of the colliding projec- 
tile and target nuclei. For other orientations of the colliding nuclei the potential is more flat 
and the depth of the pocket is smaller. Moreover, in these cases an entrance barrier and the 
minimum of the pocket of W{R) have larger absolute energies than in the case of the tip-tip 
orientation. Therefore, an excitation energy of a compound nucleus will be larger than in 
the last case. An excess of the excitation energy will decrease the survival probability of the 
evaporation residues. Thus, in the fusion of massive nuclei their mutual orientation strongly 
influences not only the capture cross section but also the probability that the compound 
nucleus survives during deexcitation. 

The existence of the window for the bombarding energy has a crucial influence on the 
fusion process. From one side a larger bombarding energy will be needed to overcome 
an intrinsic barrier {Bj^^) to form a compound nucleus. From other side an increase of 
the bombarding energy decreases the capture probability starting from some values of the 
bombarding energy because the friction force is not strong enough to provide a complete 
dissipation of the kinetic energy. 

To analyse a fusion process further, we need in a dynamic model which describes an 
evolution of a dinuclear system to compound nucleus. Below, we shall use a model developed 
in [8] . According to this model a dinuclear system evolves to compound nucleus by increasing 
its mass asymmetry. It means that driving potential (P) plays the main role in the fusion 
dynamics and a dinuclear system should overcome the Businaro-Gallone point to be fused. 
The driving potentials for the reactions which we analyse are presented in Figs. 3-5. The 
values of the barriers which should be overcome to get a compound nucleus (Bj^^) depend 
on the compound system and the reaction choice which determines the initial value of the 
mass asymmetry. These are equal to 6 MeV for ^^Ca + ^^'^Pu (Fig. 3) and 28 MeV for 
'''^'''^Ge + ^°'^Pb (Figs. 4 and 5). To overcome the barrier, a dinuclear system should have 
the corresponding excitation energy. However, the possible values of the excitation energy 
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which are defined by the amount of a dissipated energy are restricted by the framework of 
the energy window for bombarding energies leading to capture. The possible values of the 
excitation energy can be estimated and the results are shown in Fig. 6. For "^^Ca + ^^"^Pu 
reaction the excitation energy can take the values from 19 MeV up to 41 MeV which are 
larger than the barrier Bj^^ of the driving potential. In the case of '''^''''^Ge + ^°^Pb reactions, 
the excitation energy E* takes the values between 6 MeV and 16 MeV. This value is lower 
than the value of B*f^g—2S MeV for these reactions but it is larger than the quasifission 
barrier which is about (3 — 5) MeV (Fig. 2). An increase of the beam energy in order to 
obtain an adequate excitation energy does not help because dinuclear system can not be 
formed. The corresponding value of the beam energy will exceed E^ax- Thus, according 
to our calculations of a capture cross section and the model of fusion suggested in [8] , the 
compound nucleus can not be formed with a measurable cross section in the ^"^'^^Ge + ^osp^-, 
reactions. However, it is not excluded that a dinuclear system can prefer the trajectory in 
the R — r) plane for fusion different from that suggested in [8] or other mechanism of the 
compound nucleus formation like cluster transfer [11] might play an important role. 

The other question concerns the probability that the excited compound nucleus formed in 
a fusion process survives during deexcitation. An increase of an excitation energy decreases 
the influence of the shell effects on stability of a compound nucleus and decrease the fusion 
probability. However, this question is not analysed in the present paper. 

IV. CONCLUSION. 

We have analysed the partial fusion cross sections for the reactions with massive nuclei 
leading to compound nucleus with Z = 114: ^^Ca + ^''^Pu and ^^'™Ge + ^"^Pb. The main 
attention is paid to the calculations of the capture probability, which is a characteristic 
feature of an initial stage of the collision. It is shown that for the considered reactions, 
there is an energy window for the bombarding energy at which the capture cross section is 
large enough to have a physical interest. This result puts a strong limitations on the choice 
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of the bombarding energy for a given reaction. However, from other side, the excitation 
energy should be large enough to overcome an intrinsic barrier for the fusion [8] . Thus, both 
restrictions can be used to obtain an optimal choice of the projectile-target combination and 
of the bombarding energy. 
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FIGURES 

FIG. 1. The calculated capture cross section as a function of the beam energy for the ^^Ca + 
244pu (I) (full circles), ^"^Ge + ^ospb (II) (full triangles), and ^^Ge + ^ospb (HI) (open triangles) 
reactions; Bj is the Bass barrier for the reaction (i), i=I, II, and III. 

FIG. 2. The nucleus-nucleus interaction potential calculated for ^^Ge + ^'^^Pb (solid curve), 
^^Ge + ^°^Pb (dotted curve) and ^^Ca + ^^^Pu (dashed curve) reactions; Bj is the Bass barrier for 
the reaction (i), i=I, II, and III. 

FIG. 3. The driving potential for the superheavy element ^^^114. The arrow indicates an initial 
charge asymmetry which corresponds to the ^^Ca + ^^^Pu reaction. 

FIG. 4. The driving potential for the superheavy element ^^^114. The arrow indicates an initial 
charge asymmetry which corresponds to ^^Ge + ^°^Pb reaction. 

FIG. 5. The driving potential for the superheavy element ^^^114. The arrow indicates an initial 
charge asymmetry which corresponds to ''^Ge + 208p|-, reaction. 

FIG. 6. The excitation energy of a dinuclcar system formed after capture of nuclei in reactions: 
^«Ca + 244pu (I) (full circles), ^^Ge + ^^^Pb (II) (full triangles), and ^^Ge + ^^^Pb (III) (open 
triangles) as a function of the beam energy in the center of mass system. 
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